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<N ■ Abstract 

The quantum theory of spin light (electromagnetic radiation emitted by a Dirac massive 
neutrino propagating in dense matter due to the weak interaction of a neutrino with back- 
ground fermions) is developed. In contrast to the Cherenkov radiation, this effect does not 
disappear even if the medium refractive index is assumed to be equal to unity. The formulas 
for the transition rate and the total radiation power are obtained. It is found out that radiation 



o 



of photons is possible only when the sign of the particle helicity is opposite to that of the 
q \ effective potential describing the interaction of a neutrino (antineutrino) with the background 

medium. Due to the radiative self-polarization the radiating particle can change its helicity. 
As a result, the active left-handed polarized neutrino (right-handed polarized antineutrino) 
Qh| converting to the state with inverse helicity can become practically "sterile". Since the sign of 

Q_i| the effective potential depends on the neutrino flavor and the matter structure, the spin light 

Jjn \ can change a ratio of active neutrinos of different flavors. In the ultra relativistic approach, 

the radiated photons averaged energy is equal to one third of the initial neutrino energy, and 

• rH . 

^ , two thirds of the energy are carried out by the final "sterile" neutrinos. 

a' 

A Dirac massive neutrino has non-trivial electromagnetic properties. In particular, it pos- 
sesses non-zero magnetic moment [1]. Therefore a Dirac massive neutrino propagating in dense 
matter can emit electromagnetic radiation due to the weak interaction of a neutrino with back- 
ground fermions [2, 3]. As a result of the radiation, neutrino can change its helicity due to the 
radiative self-polarization. In contrast to the Cherenkov radiation, this effect does not disappear 
even if the refractive index of the medium is assumed to be equal to unity. This conclusion is 
valid for any model of neutrino interactions breaking spatial parity. The phenomenon was called 
the neutrino spin light in analogy with the effect, related with the synchrotron radiation power 
depending on the electron spin orientation (see [4]). 
*E-mail: lobanov@phys.msu.ru 
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The properties of spin light were investigated basing upon the quasi-classical theory of radia- 
tion and self-polarization of neutral particles [5,6] with the use of the Bargmann-Michel-Telegdi 
(BMT) equation [7] and its generalizations [8,9]. This theory is valid when the radiated photon 
energy is small as compared with the neutrino energy, and this narrows the range of astrophysical 
applications of the obtained formulas. 

In the present paper, the properties of spin light are investigated basing upon the consistent 
quantum theory, and this allows the neutrino recoil in the act of radiation to be considered for. 
The above mentioned restriction is eliminated in this way. 

On the other hand, the detailed analysis of the results of our investigations shows that the 
features of the effect depend on the neutrino flavor, helicity and the matter structure [10]. This 
fact leads to the conclusion that the spin light can initiate transformation of a neutrino from the 
active state to a practically "sterile" state, and the inverse process is also possible. 

When the interaction of a neutrino with the background fermions is considered to be coherent, 
the propagation of a massive neutrino in the matter is described by the Dirac equation with 
the effective potential [11, 12]. In what follows, we restrict our consideration to the case of a 
homogeneous and isotropic medium. Then in the frameworks of the minimally extended standard 
model, the form of this equation is uniquely determined by the assumptions similar to those 
adopted in [13]: 

id- l -f{l+^)-m^W v = Q. (1) 

The function f M is a linear combination of fermion currents and polarizations. The quantities 
with hats denote scalar products of Dirac matrices with 4- vectors, i.e. , a = y^a^. 

If the medium is at rest and unpolarized then f = 0. The component f° calculated in the first 
order of the perturbation theory is as follows [14-16]: 

f = V2G F { ^ [hv + - 2Q {f) sin 2 9 W ) (n f - n f )\. (2) 

Here, n f , rif are the number densities of background fermions and anti-fermions, Q ( f) is the elec- 
tric charge of the fermion and T^p is the third component of the weak isospin for the left-chiral 
projection of it. The parameter I ev is equal to unity for the interaction of electron neutrino with 
electrons. In other cases I ev = 0. Summation is performed over all fermions / of the background. 

Let us obtain a solution of equation ([T]). Since function f M = const, equation © commutes 
with operators of canonical momentum id^. However, the commonly adopted choice of eigenval- 
ues of this operator as quantum numbers in this problem is not satisfactory. Kinetic momentum 
components of a particle, related to its group 4-velocity u M by the relation q u = niyU^, q 2 = m 2 ., 
are more suitable to play the role of its quantum numbers. This choice can be justified, since it is 
the particle kinetic momentum that can be really observed. 
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The explicit form of the kinetic momentum operator for the particle with spin is not known 
beforehand, and hence, in order to find the appropriate solutions, we have to use the correspon- 
dence principle. 

It was shown in [8] that when the effects of the neutrino weak interaction are taken into ac- 
count, the Lorentz invariant generalization of the BMT equation for spin vector S M is as follows: 

S" = 2/u Q { (F» v + 0' V )S V - u»u v (F vA + G VA )S A ), (3) 

where 

cr v = ^ vpA f P u A , (4) 

and a dot denotes the differentiation with respect to the proper time r. 

Let us introduce the quasi-classical spin wave functions. Such wave functions can be con- 
structed as follows [6,9]. Suppose the Lorentz equation is solved, i. e. the dependence of particle 
coordinates on proper time is found. Then the BMT equation transforms to ordinary differen- 
tial equation, whose resolvent determines a one-parametric subgroup of the Lorentz group. The 
quasi-classical spin wave function is represented by a spin-tensor, whose evolution is determined 
by the same one-parametric subgroup. 

In the case when the effect of an external electromagnetic field can be neglected as com- 
pared with the effect of the neutrino interaction with the background matter, the equation for the 
neutrino quasi-classical wave function ¥(r) is 

W(t) = i/ioy^G^W Y(t), (5) 

where *G MV = - l l 2 e^ vpA G pA is a tensor dual to G MV . Obviously, the quasi-classical density matrix 
of a polarized neutrino takes the form 

p(T,r') = X - U(T,T )(q(To) + m)(l- 7 5 S(T ))U- l (T',ToX (6) 
where U(t, t ) is the resolvent of equation ©, and the equation for U(r, To) is 

U(t, t ) = l - y 5 (fu - uf) U(r, r ). (7) 

We note that the operator U(t, t ) is defined up to a phase factor e~ lF<x \ with the derivative of the 
exponent with respect to the proper time is equal to zero: 

F{x) = 0. (8) 

Let us choose the solution of equation ([Hi in the form [9] 

W(x) = U(r(x))W (x), (9) 
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where ¥ is a solution of the Dirac equation for a free particle 

W (x) = e~ i{qx \q + m v )(l - y 5 S ¥°. 
Here is constant bispinor and ¥ (x) normalized by the condition 

%(x)W (x) = 2my. 
Substitution of the expression © in eq. © results in the relation 

{§ + (OF) - X -f + l -y 5 f + i y 5 N (fu - uf) - m y | e~ iF(x) U (r(x)) ^ = 0, 



(10) 



(11) 



where N M = c^t. Since the commutator [q, U] = 0, and the matrix U is nondegenerate, then for 
this relation to hold the following condition is required 



(3F) - -/ + - 7 5 f + - 7 5 n (fu - uf) = o. 



r i' j 4 

It is easy to find out that the equation (fT21) is valid only if 

d A F = l -f\ e» vp %f v u p = 0, (1 - (Nu)) f + (Nf)u A = 0. 



From two latter equations it follows that 



nfu) - u"f 2 

(fuf-fu 2 ■ 



So f M = const, then 



and we can write 



Here 



T = (Nx), F = -(fx), 

U(x) = e- iifx)/1 e^A^. 
<r=±i 

A ? = -[l-£y 5 S tp q/m v \, £±1, 
are spin projection operators with eigenvalues £ ± 1 respectively, and 

t 1/ r \? 7T~i WW-/ 2 ^) c , q»(fq)/m v -Pm v 
<p = -^m-Prn v = 2 ^____, S tp = V (/g)2 _ /2m , 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 
(18) 



From the obtained formulas it follows that the eigenvalues of the operator of canonical mo- 
mentum ie^ are 



1 + 



^f 2 



2 V(/4) 2 " 



+ 



1 - 



^(fqf-fhnl) 



(19) 



The dispersion law follows from eq. (fT9"l) in the form 

P 2 = m 2 v + (Pf) -f 2 /2-£ ^j((Pf)-f2/2f-p m 2. 



(20) 



If the medium is at rest and unpolarized then the neutrino total energy and canonical momen- 
tum are determined by the formulas 

e = q° + f/2, P = qA 9 f, (21) 

where A^ = 1 + £/°/2|q|, and 

S1 P = — (iqUVlql), (22) 

i. e. the eigenvalues £ = +1 determine the helicity of the particle. Consequently, the dispersion 
law is 

s = yl(A\P\-£f>/2f+ml + f/2, (23) 
where A = sign(A^). Obviously 

de q 
dP = ^ 

is the particle group velocity. 

The relation (1231 differs those used in previous papers (see, for example, [18]) by the mul- 
tiplier A. This is due to the fact that, in these papers the projection of the particle spin on the 
canonical momentum P and not the helicity of the particle was used as the spin quantum number 
The helicity is the projection of the spin on the direction of its kinetic momentum [19-21], 
because the rest frame of the particle is determined by the condition that its group velocity is 
equal to zero. In our problem the directions of canonical and kinetic momenta, generally speak- 
ing, are different, and hence, the projection of particle spin on the canonical momentum does not 
coincide with its helicity. 

From formulas d2TT) . it is seen that if the energy is expressed in terms of the kinetic momen- 
tum, then it does not depend on the particle helicity, while the particle canonical momentum is a 
function of the helicity. Therefore, the statement of the authors of [17], i.e., that the radiation of 
photons in the process of the spin light emission takes place due to neutrino transitions from the 
"exited" helicity state to the low-lying helicity state in matter, is not correct. 

Let us consider the process of emitting photons by a massive neutrino in unpolarized matter 
at rest. In this case, the orfhonormalized system of solutions for equation © is: 

¥{x) = -M e -<w o /2>.vy^ + my)(1 _ tf§ ¥ o (24) 
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The formula for the spontaneous radiation transition probability of a neutral fermion with 
anomalous magnetic moment fi is 1 : 

P = --i- fd 4 xd 4 y f ^±S(k 2 )S(q 2 -m 2 v )X 

2p°J "J (2n) 6 w vJ (25) 

xSpir^giix, y; p, £i)r v (y)g f (y, x; q, 0}g^ h (x, y; k). 

Here, g t {x, y; p), Q/(y, x; q) are density matrices of the initial (z) and final (/) states of the fermion, 
(fj h (x, y; k) is the radiated photon density matrix, P 1 = - ^[An^cr^ky is the vertex function. The 
density matrix of longitudinally polarized neutrino in the unpolarized matter at rest constructed 
with the use of the solutions of equation (|TJ has the form 

g(x,y; p, = \^{p + m v )(l - tfS^W+J*!™**^. (26) 

After summing over photon polarizations 2 and integrating with respect to coordinates we 
obtain the expression for the transition rate under investigation: 

W = ^J d -^~ W 2 )8(q 2 ~ m 2 v )6(p° - q° - V(pA K , - q\ (f - k)T(p, q), (27) 

where 

T(p, q) = 4A 2 p( .A 2 q(f {(pk)(qk) - [/c°|p| - p°(pk)/|p|] [*°|q| - <Aqk)/|q|]) . (28) 
After integrating over k, k°, |q| we obtain the spectral-angular distribution of the final neutrino 

2 P ° 

W = ~^ f ^p\\f dq ° A '^' f d0X 

m v 

Xd((p° - q ) 2 + 2|p||q|A K A^ cos & v - |p| 2 A^. - |q| 2 A^) x (29) 



where 



x{(/°/2) 2 [^|p||q| + ami - A )] 2 +[(/ /2)(^°|p| - <T//lql) + m 2 v (p° - q )] 

Iql = ^(q°) 2 -m 2 v . 

It is convenient to express the results of integrating over angular variables using dimension- 
less quantities. Introducing the notations 

x = q Q /m v , y = p°/m v , d = \f\/2m v , fa = fa sign(/°) (30) 



'in the expression for the radiation energy <5, the additional factor k, i.e. the energy of radiated photon, appears 
in the integrand. 

2 We do not consider the polarization of spin light photons here. In the quasi-classical approximation, this problem 
was investigated in [17]. 
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we have 



2 3 

w 



if l = h 



xed) ye[l,y ), 

xe[a)i,a) 2 ] ye[y ,yi), 

xg[1,oj 2 ] ye [71,72), 

x e ye [72, 00), 



if £ = = -1, and 



if £• = -1,?/ = 1. 

Here 
where 



and 



Zi = 7 + Vt 2 _ 1 _ 2d, 
Z2 = 7 - Vt 2 _ 1 + 2 <i, 

y = Vl+J 2 , 

ri = i{(l+2J) + (l+2dr 1 }, 

72 = ^j(l-2d) + (l-2d)- 1 } d<l/2, 



(3D 



+ [ 7 -x + d (fa Vr 2 -i - r/r Vx 2 - r 

The integration bounds in the formula (OTT) are 

x£0 ye [l,oo), (32) 



(33) 



xe0 ye [1,70, 

xe[l,6>!] ye[y!,y 2 ), (34) 
x e [aj 2 ,&>i] 7 e [72, 00). 



iO\ = i(zi + V)> C0 2 = ^(z 2 +22 1 )' ( 35 ) 



(36) 



(37) 



72 = 00 d > 1/2. 

The integration is carried out elementary. The transition rate under investigation is defined as 



2 3 

% = "V^ (1 + 0) D©(r - n) + z fe, -D0(r - 72)] + 
+(l - [Z(zi, l)0(7i - 7) + Zfe, -1)0(72 - y)] 0(y - yo))d - ft). 
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(38) 



Here 

Z fc ?/) = ^TTT) i lnz b 2 + d + d 2 + l/2} + 

+\{z 2 ~ z~ 2 ) [d 2 (2y 2 -\)+d Vr 2 - 1 + l/2] + 

+ ^(z-z') 2 [2djy^l + l]dy- (39) 
-(z-z~ 1 )[d 2 + dyjy 2 - 1 + l]y- 

-^(z + z-'-^^Vr^ + r 2 ]^}- 

Therefore, the transition rate after summation over polarizations of the final neutrino becomes 



2 3 

w lf=i + w^ = _! = ^(i - £){z(zi, l) + z(z 2 , -i))0( r - ro ). (40) 



2 

If dy <sc 1, then expression (l38l) leads to the formula 



% = — ^— (r z - D J/z (i - + (41) 



,2 _ 1 ^3/2 



obtained in the quasi-classical approximation in [3]. 

In the ultra-relativistic limit (7 » 1, dy » 1), the transition rate is given by the expression 

W lf =^mld 2 y{\-ld{l+l f ). (42) 
Let us consider now the radiation power. If we introduce the function 

Z{z,lf) = yZ{z,lf)-Y{z,lf), (43) 



where 



y fc & = ^rij I " ln z K + d + 1] r - 

-5 (z 2 - z" 2 ) [d 2 + d Vr 2 - 1 + 1] r+ 
+^ (z - z' 1 ) 3 [d 2 (2y 2 -i) + d Vr 2 -i + 1/2] + 

+- (z - z" 1 ) [2J 2 y 2 + 2d V7 2 - 1 + 7 2 + l] + 
12 



(44) 



+^ ((z + z" 1 ) 3 - s) [2d Vr 2 "^" + 1] ^7- 



then the formula for the total radiation power can be obtained from (I3£b . (EfUb by the substitution 
Z(z, If) — » Z(z, if). It can be verified that if dy <sc 1 then the radiation power is 



32u 2 m 4 d 4 

h } = 3 V (/ - W ~ 6X1 + fr). (45) 



.2 i \2 



This result was obtained in the quasi-classical approximation in [2]. In the ultra-relativistic limit, 
the radiation power is equal to 

1 



/^-^y r 2 (i-( i )(i+( / ). (46) 

It can be seen from equations d42l and d46T) that in the ultra-relativistic limit the averaged energy 
of emitted photons is (s y ) = s v /3. It should be pointed out that the obtained formulas are valid 
both for a neutrino and for an antineutrino. The charge conjugation operation leads to the change 
of the sign of the effective potential and the replacement of the left-hand projector by the right- 
hand one in the equation O. Thus the sign in front of the y 5 matrix remains invariant. 

Using eq. d27t . it is possible to find the dependence of the radiated photon energy on the 
angle § y between the direction of the neutrino propagation and the photon wave vector: 

* = 2d PX - dh . (47) 
m v (X + d/y)(X-d/y) V J 



Here fi = -\jy 2 - l/y is the neutrino velocity and X = 1 — (J3 — d/y) cos # r . In the quasi-classical 
approximation, this formula reduces to the relation 

- = r^r- < 48 > 

m v 1 -pcosVy 

which follows from the results of [3] after Lorentz transformation to the laboratory frame. 

The following conclusions can be made from the obtained results. A neutrino (antineutrino) 
can emit photons due to coherent interaction with matter only when its helicity has the sign 
opposite to the sign of the effective potential f°. Otherwise, radiation transitions are impossible. 
In the case of low energies of the initial neutrino, only radiation without spin-flip is possible and 
the probability of the process is very low. At high energies, the main contribution to radiation 
is given by the transitions with the spin-flip, the transitions without spin-flip are either absent or 
their probability is negligible. This leads to total self-polarization, i. e. the initially left-handed 
polarized neutrino (right-handed polarized antineutrino) is transformed to practically "sterile" 
right-handed polarized neutrino (left-handed polarized antineutrino). For "sterile" particles, the 
situation is opposite. They can be converted to the active form in the medium "transparent" for 
the active neutrino. 
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With the use of the effective potential calculated in the first order of the perturbation theory 
©, the following conclusions can be made. If the matter consists only of electrons then, in the 
framework of the minimally extended standard model in the ultra-relativistic limit (here we use 
gaussian units), we have for the transition rate 



and for the total radiation power 



96 h\iu B ) \m e c 



k f = ^\^\ l^fl (1-6X1 (50) 



Here s v is the neutrino energy, /i B = e/2m e is the Bohr magneton, or is the fine structure constant, 
m e is the electron mass and Gp = Gp(l + 4 sin 2 # w ), where Gp, #w are the Fermi constant and the 
Weinberg angle respectively. Thus, after the radiative transition, two thirds of the initial active 
neutrino energy are carried away by the final "sterile" one. 

At the same time, as it can be seen from eq. ©, a muon neutrino in the electron medium 
does not emit any radiation. Moreover, a muon neutrino does not emit radiation in an electrically 
neutral medium, when the number density of protons is equal to the electron number density. 
And an electron neutrino can emit radiation if the electron number density is greater than the 
neutron number density. An example of such medium is provided by the Sun. Therefore the spin 
light can change the ratio of active neutrino of different flavors. 

It is obviously that the above conclusions change to opposite if the matter consists of an- 
tiparticles. Therefore the neutrino spin light can serve as a tool for determination of the type of 
astrophysical objects, since neutrino radiative transitions in dense matter can result in radiation of 
photons of super-high energies, even exceeding the GZK cutoff. Indeed, the neutron medium is 
"transparent" for all active neutrinos, but an active antineutrino emits radiation in such a medium, 
the transition rate and the total radiation power can be obtained from equations d4*9t and (l50l) af- 
ter substitution G F — » G F ,n e — > n n . If the density of the neutron star is assume to be n « 10 38 , 
the transition rate is estimated as 

/ ^ 2 

22 S v I 



W= 10 z/ ^ — , (51) 
£gzk \I^b I 

where s GZ k = 5 x 10 19 eV is GZK cutoff energy. Although the transition rate determined by eq. 
(BTT) is extremely low, this effect can still serve as one of a possible explanations of the cosmic 
ray paradox. 

The spin light can also be important for the understanding of the "dark matter" formation 
mechanism in the early stages of evolution of the Universe. 
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When the present paper was already submitted for publication, we came across an article [22], 
where the spin light theory was also considered. The formulas of [22] in the ultra-relativistic limit 
of physical interest reproduce the results for the transition rate and the total power of spin light 
already obtained in our earlier publication [10]. 

The author is grateful to V.G. Bagrov, A.V. Borisov, and V.Ch. Zhukovsky for fruitful discus- 
sions. 
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References 

[1] K. Fujikawa, R. Shrock, Phys. Rev. Lett. 45 (1980) 963. 

[2] A. Lobanov, A. Studenikin, Phys. Lett. B 564 (2003) 27, hep -ph/02123"93l 

[3] A. Lobanov, A. Studenikin, Phys. Lett. B 601 (2004) 171, |astKv ph/0408026. 



[4] V.A. Bordovitsyn, I.M. Ternov, V.G. Bagrov, Physics - Uspekhi 38 (1995) 1037, (Uspekhi 
Fiz. Nauk 165 (1995) 1084). 



[5] A.E. Lobanov, O.S. Pavlova, Phys. Lett. A 275 (2000) 1, hep-ph/0006291 



[6] A.E. Lobanov, hep-ph/03 1 1021 



[7] V. Bargmann, L. Michel, V. Telegdi, Phys. Rev. Lett. 2 (1959) 435. 

[8] A. Lobanov, A. Studenikin, Phys. Lett. B 515 (2001) 94, |hep-ph/0106101| 

[9] A.E. Lobanov, O.S. Pavlova, Moscow Univ. Phys. Bull. 54(4) (1999) 1, (Vestn. MGU. Fiz. 
Astr. 40(4) (1999) 3). 



[10] A.E. Lobanov, hep-ph/041 1342 



[11] L. Wolfenstein, Phys. Rev. D 17 (1978) 2369. 

[12] S.P Mikheyev, A.Yu. Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913, (Yad. Fiz. 42 (1985) 
1441). 

[13] L.L. Foldy, Phys. Rev. 87 (1952) 688. 



11 



[14] P.B. Pal, T.N. Pham, Phys. Rev. D 40 (1989) 259. 

[15] J.F. Nieves, Phys. Rev. D 40 (1989) 866. 

[16] D. Notzold, G. Raffelt, Nucl. Phys. B 307 (1988) 924. 



[17] A. Studenikin, A. Ternov, hep-ph/04 10297, hep-ph/04 12408 



[18] J. Pantaleone, Phys. Lett. B 268 (1991) 227. 

[19] A. A. Sokolov, I.M. Ternov, Synchrotron radiation, Berlin, Academic Verlag (1968). 

[20] V.G. Bagrov, D.M. Gitman, Exact Solutions of Relativistic Wave Equations, Dordrecht, 
Kluwer (1990). 

[21] I.M. Ternov, Introduction to Physics of Spin of Relativistic Particles [in Russian], Moscow, 
Moscow State Univ. (1997). 



[22] A. Grigoriev, A. Studenikin, A. Ternov, hep-ph/050223 1 



12 



